Mon. Not. R. Astron. Soc. 000, 1-14 (2009) Printed 24 March 2009 (MN KTbX style file v2.2) 



Asteroseismic analysis of the roAp star a Circini: 84 days of 
high-precision photometry from the WIRE satellite"^ 

H. Bruntt\ D. W. Kurtz^, M. S. Cunha^ 1. M. Brandao^'^, G. Handler^ T. R. Bedding^ 
T. Medupe^'^, D. L. Buzasi^ D. Mashigo^ 1. Zhang^ F. van WyF 

^Sydney Institute for Astronomy, School of Physics A29, University of Sydney, 2006 NSW, Australia 

^Jeremiah Horrocks Institute of Astrophysics and Supercomputing, University of Central Lancashire, Preston PRI 2HE, UK 

^Centro de Astrofisica da Universidade do Porto, Rua das Estrelas, 4150-Porto, Portugal 

'^Departamento de Matemdtica Aplicada, Faculdade de Ciencias, Universidade do Porto, 4169 Porto, Portugal 

^Institut fur Astronomic, Universitdt Wien, TUrkenschanzstrasse 17, A- 1180 Wien, Austria 

^Department of Astronomy, University of Cape Town, Rondebosch 7700, South Africa 

^ South African Astronomical Obsen>atory, P.O. Box 9, Obser\>atory 7935, South Africa 

^ Department of Astronomy, University of Washington, Seattle WA 98195, USA 



Accepted XXX; Received YYY; in original form ZZZ 



ABSTRACT 

We present a detailed study of the pulsation of a Circini, the brightest of the rapidly oscillating 
Ap stars. We have obtained 84 days of high-precision photometry from four runs with the 
star tracker on the WIRE satellite. Simultaneously, we collected ground-based Johnson B 
observations on 16 nights at the South African Astronomical Observatory. In addition to the 
dominant oscillation mode at 2442 /iHz, we detect two new modes that lie symmetrically 
around the principal mode to form a triplet. The average separation between these modes is 
A/ = 30.173 ± 0.004 /iHz and they are nearly equidistant with the separations differing by 
only 3.9 nHz. We compare the observed frequencies with theoretical pulsation models based 
on constraints from the recently determined interferometric radius and effective temperature, 
and the recently updated Hipparcos parallax. We show that the theoretical large separations 
for models of a Cir with global parameters within the 1 a observational uncertainties vary 
between 59 and 65 /iHz. This is consistent with the large separation being twice the observed 
value of A/, indicating that the three main modes are of alternating even and odd degrees. 
The frequency differences in the triplet are significantly smaller than those predicted from our 
models, for all possible combinations of mode degrees, and may indicate that the effects of 
magnetic perturbations need to be taken into account. The WIRE light curves are modulated 
by a double wave with a period of 4.479 d, and a peak-to-peak amplitude of 4mmag. This 
variation is due to the rotation of the star and is a new discovery, made possible by the high 
precision of the WIRE photometry. The rotational modulation confirms an earlier indirect 
determination of the rotation period. The main pulsation mode at 2442 /iHz has two sidelobes 
split by exactly the rotation frequency. From the amplitude ratio of the sidelobes to the central 
peak we show that the principal mode is consistent with an oblique axisymmetric dipole mode 
{I = 1,TO = 0), or with a magnetically distorted mode of higher degree with a dominant 
dipolar component. 

Key words: stars: individual: a Cir - stars: oscillations - stars: rotation - stars: fundamental 
parameters - stars: variables: other 



1 INTRODUCTION 

a Circini (HR5463, HD 128898, HIP71908; V = 3.2) is the 
brightest of the rapidly oscillating chemically peculiar A-type 
(roAp) stars (see reviews by Kurtz & Martinez 2000 and Kurtz 

* Based on observations at the South African Astronomical Observatory. 



et al. 2004). It was among the first of the class to be discov- 
ered (Kurtz 1982), showing photometric variations with a period 
Pp\A = 6.8 min and an amplitude of a few mmag in Johnson B 
(Kurtz &. Cropper 1981). The multi-site photometric campaign by 
Kurtz et al. (1994) revealed four low amplitude modes in addition 
to the main mode. The roAp stars pulsate in high overtones and 
are described by the oblique pulsator model (Kurtz 1982; Shiba- 
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Table 1. List of photometric observations from the WIRE satellite and the 
0.5-m and 0.75-m telescopes at SAAO. In the last column we list the point- 
to-point scatter (cr) in the time series. 



Source 


Observing dates 


'^obs 


Data 








[days] 


points 


[mmag] 


WIRE 


2000 Aug. 18-Sep. 29 


42.2 


39479 


0.93 




2005 Feb. 16-25 


8.3 


10 857 


1.04 




2006 Feb. 10-Mar. 15 


31.0 


29616 


0.91 




2006 July 7-31 


24.3 


27 273 


0.61 


SAAO 


2006 Feb. 7-Feb. 19 


12.0 


1 802 


2.4 




2006 June 23-July 24 


31.2 


2 330 


2.0 



hashi & Takata 1993; Takata & Shibahashi 1995; Bigot & Dziem- 
bowski 2002; Saio 2005). In this model the magnetic axis is in- 
clined to the rotation axis. The oscillation modes follow the mag- 
netic field, therefore the deformation of the star is observed from 
different perspectives. This was observed in a Cir by Kurtz et al. 
(1994), who found a significant modulation of the amplitude of the 
main mode, making it possible to determine the rotation period to 
be P,ot = 4.4790 ± 0.0001 d. 

Simultaneous excitation of several high-order modes in roAp 
stars makes them ideal targets for detailed asteroseismic studies. 
Asymptotic theory predicts a regular pattern of frequencies, domi- 
nated by the large separation, which depends mainly on the mean 
density of the star. Asteroseismology of roAp stars makes it pos- 
sible to probe their internal properties, including aspects of the 
structure and magnetic fields (Cunha et al. 2003). Previous detailed 
studies of multi-mode roAp stars have in some cases been lim- 
ited because the fundamental parameters are poorly known (e.g. 
Cunha et al. 2007). For example, the effective temperature deter- 
mined from different photometric indices or spectroscopic analy- 
ses typically give results spanning 500 K or more (Netopil et al. 
2008). This large range may be explained, in part, by the chemi- 
cal peculiarities and strong magnetic fields in roAp stars, resulting 
in heavy line blanketing in some regions of the spectra. Matthews 
et al. (1999) presented observational evidence of systematic errors 
in the effective temperatures of a dozen roAp stars, based on the 
measurements of their large separations. 

In our companion paper (Bruntt et al. 2008, Paper I), we used 
interferometric measurements with the Sydney University Stellar 
Interferometer (SUSI) to determine the radius of a Cir. Combined 
with a determination of the bolometric flux from calibrated spectra, 
we found the temperature to be 7420 ± 170 K. Unlike previous de- 
terminations of Teff, which range from 7670 to 8440 K, our deter- 
mination is nearly independent of atmospheric models - only the 
limb-darkening correction depends on the adopted grid of model 
atmospheres, and weakly so. 

Here we report observations of aCir, made simultaneously 
with the WIRE satellite and from the ground. The high-precision 
satellite photometry has allowed us to detect the large separation 
in Q Cir for the first time. This new result, combined with the con- 
straints on the radius and T^s from Paper I, allow us to make an 
asteroseismic investigation based on theoretical models. 

2 OBSERVATIONS 

Soon after the launch of the WIRE satellite on 4 March 1999, the 
hydrogen coolant for the main infrared camera was lost. However, 



Table 2. List of photometric observations from the 0.5-m and 0.75-m tele- 
scopes at SAAO. The first and second columns give the initials of the ob- 
server and the telescope used. The duration in the fourth column is from 
start to end of observations for each night which sometimes includes gaps 
due to poor weather. The fifth column gives the number of averaged 40-s 
integrations obtained and the last column is the standard deviation of one 
observation with respect to the mean for the night. 



Obs. 


Tel 


Date 


-^obs 


Data 


cr 






2006 


[h] 


points 


[mmag] 


GH 


0.5-m 


07-08 Feb. 


0.58 


49 


3.5 


GH 


0.5-m 


12-13 Feb. 


0.95 


83 


3.4 


GH 


0.75-m 


13-14 Feb. 


3.04 


258 


1.8 


GH 


0.5-m 


13-14 Feb. 


0.45 


40 


2.0 


GH 


0.75-m 


14-15 Feb. 


3.97 


341 


2.4 


GH 


0.75-m 


15-16 Feb. 


2.58 


224 


1.8 


GH 


0.75-m 


16-17 Feb. 


3.95 


264 


2.3 


GH 


0.75-m 


17-18 Feb. 


3.49 


302 


2.7 


GH 


0.75-m 


19-20 Feb. 


2.88 


241 


2.5 


total 






21.89 


1802 




DM 


0.5-m 


23-24 Jun. 


5.45 


381 


2.0 


DM 


0.5-m 


24-25 Jun. 


5.45 


445 


1.7 


DM 


0.5-m 


07-08 Jul. 


5.22 


418 


1.9 


DM 


0.5-m 


11-12 Jul 


3.89 


88 


2.1 


DM 


0.5-m 


17-18 Jul. 


4.79 


396 


1.8 


DM 


0.5-m 


18-19 Jul. 


4.79 


389 


2.0 


IZ 


0.5-m 


24-25 Jul. 


4.16 


213 


1.9 


total 






33.75 


2330 





the star tracker was used to observe around 240 stars with apparent 
magnitudes V < Q until communication with the satellite failed 
on 2006 October 23 (Buzasi 2002; Bruntt 2007; Bruntt & South- 
worth 2008). The 52 mm aperture illuminated a SITe CCD cam- 
era and windows centered on the five brightest stars in the field 
were acquired at a cadence of 2 Hz. Each WIRE field was typically 
monitored for 2 to 4 weeks and some targets, including a Cir, were 
observed in more than one season. WIRE switched between two 
targets during each orbit in order to best avoid the illuminated face 
of the Earth. Consequently, the duty cycle was typically around 20 
to 40 percent, or about 19 to 38 minutes of continuous observa- 
tions per orbit. The time baseline of the observations analysed here 
is almost six years, during which the orbital period of the satellite 
decreased from 96.0 to 93.3 minutes. 

The roAp star a Cir was observed during four runs with WIRE. 
The first data set in 2000 September was not optimal in that several 
long gaps are present and the sub-pixel position on the CCD was 
less stable than the succeeding runs. The second run in 2005 Febru- 
ary lasted 8 days, and was followed by runs of 31 and 24 days in 
2006 February and July. A log of the observations is given in Ta- 
ble 1. Excluding the gaps in the time series, a Cir was observed for 
84 days, making it one of the most observed stars with WIRE. 

To be able to compare the WIRE results with previous ground- 
based studies, we collected additional ground-based photometric 
time-series. During both runs with WIRE in 2006 we made simul- 
taneous photometric observations in the Johnson B filter from the 
South African Astronomical Observatory (SAAO). We used the 
Modular Photometer on the 0.5-m telescope and the University of 
Cape Town Photometer on the 0.75-m telescope. All observations 
were made through a Johnson B filter, plus a neutral density filter 
of 2.3 mag on the 0.5-m telescope and 5 mag on the 0.75-m tele- 
scope. The neutral density filters were needed to keep the count 
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Figure 1. The four light curves of a Cir from WIRE and the two light curves 
from SAAO. The different zero points in time are indicated in each panel. 
The bottom panel shows an expanded view of three days where simultane- 
ous observations with WIRE and SAAO were made. 



rates for this very bright star at levels that would not damage the 
photomultiplier tubes, or result in large dead-time corrections. The 
resultant count rates were about 380 000 s"^ for the 0.5-m tele- 
scope and 150000s~^ for the 0.75-m telescope. A detailed log of 
the observations from SAAO is given in Table 2. 

The light curves from WIRE and SAAO are shown in Fig. 1. 
The bottom panel shows an expanded view of simultaneous obser- 
vations. In each panel we give the time zero point to as the Helio- 
centric Julian Date (HJD). 

2.1 Data reduction 

The WIRE data set consists of 3.3 million windows extracted by 
the on-board computer from the 512x512 pixels CCD star tracker 
camera. Each window is 8 x 8 pixels and aperture photometry was 
carried out using the pipeline described by Bruntt et al. (2005). 
The data collected within 15 seconds are binned to decrease the 
number of data points to 107 225 (cf. Table 1). We note that a Cir 
is a visual binary with a separation of 15"6. Since the pixel size 
of the WIRE CCD is 60", the B component is entirely within the 
applied photometric aperture. As noted by Kurtz et al. (1994) the 
contamination is negligible, since the difference in magnitude is 
AV — 5.05. With the Stromgren indices of the B-component from 



Sinachopoulos (1989) and using the TEMPLOGG software (Rogers 
1995; Kupka & Bruntt 2001) for the calibration, we find Tcff ~ 
4660 ± 200 K, log 3 = 4.6 ± 0.2, and [Fe/H] = -0.26 ± 0.10, 
suggesting that it is a K5 V star. In the following, we use a Cir to 
mean the A component of the system. 

The raw light curves contain a small fraction (1 to 10%) of 
outlying data points. These points were identified and removed 
based on the measured position, background level, and width of 
the stellar image (PSF). The light curves from Sep'OO, Feb'05 and 
Feb' 06 have a high amount of scattered light towards the end of 
each orbit, increasing the background level from about 20 to 1400 
ADU. In these light curves we found a significant correlation be- 
tween the background and the measured flux, which we removed 
by a spline fit. The light curve from Jur06 has low background, 
and in this case we found a correlation of the measured flux and 
the width of the PSF, which was removed by subtracting a second- 
order polynomial. 

The SAAO data were obtained in high-speed mode using con- 
tinuous 10-s integrations. With continuous observations of the tar- 
get star, no comparison stars were observed. Sky observations were 
made as needed, based on conditions throughout the observing run. 
The data were corrected for dead-time losses, sky subtracted and 
the mean extinction was removed. They were then averaged four 
points at a time into 40-s integrations. Some remaining low fre- 
quency variability, caused by transparency variations, was removed 
with high-pass filtering that did not affect the pulsation frequen- 
cies of interest in the data. This latter step results in white noise 
across the frequency spectrum, and hence better estimates of errors 
on amplitudes and phases than would be the case if low frequency 
variance were left in the data. The last column of Table 2 gives the 
standard deviation for each night of observations. The contribution 
to the standard deviation from the actual pulsations is small, so it a 
good measure of the quality of the night. On some nights, weather 
interruptions resulted in gaps in the data, which is reflected in the 
number of data points and duration of each night in Table 2. 



3 LIGHT CURVE ANALYSIS 

The light curves of a Cir show variations on two different time 
scales. The double wave seen in the WIRE data in Fig. 1 has a pe- 
riod of 4.479 d and is due to rotation, seen here for the first time, 
while the periods of the rapid oscillations lie in the range 6.5 to 
7.0 min, or three orders of magnitude faster. In the following Sec- 
tions we will analyse these variations in detail. 



3.1 Rotational period and inclination angle 

To measure the rotation period, we fitted the observed double 
wave with two sinusoids with frequencies /a and fs ~ 2 x /a, 
yielding a period of P^ot = 4.4792 ± 0.0004 d. In this analysis 
we first subtracted the high-frequency oscillations (see Sect. 3.2). 
In Fig. 2 we show the binned lights curves from Feb'06 and 
Jur06, phased with the rotation period with a zero point in time 
trot ~ HJD245 3937.2086. The peak-to-peak variation is only 
about 4 mmag and it is the excellent precision and low frequency 
stability of the WIRE data that have allowed us to make the first 
direct measurement of the rotation of a Cir. Ground-based obser- 
vations have not been able to do this, partly because of the low 
amplitude of the variations and partly because of the challenge of 
making differential photometric observations on such a bright star. 
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Figure 2. Binned light curves of a Cir from WIRE phased with the rotation 
period. Triangle symbols are data from Feb'06 and circles are from Jur06. 
Vertical lines are shown to guide the eye. The zero point in time for the 
phase is i^ot = HJD245 3937.2086. 

From the rotation period and the interferometrically deter- 
mined radius of i?/RQ = 1.967 ± 0.066 (Paper I), we find the 
rotational velocity to be v^ot — 21.4 ± 0.7 km . We can thus es- 
timate the inclination of the rotation axis from the measured v sin i. 
The average value from Kupka et al. (1996) and Balona & Laney 
(2003) is 11 sin i — 13.0 ± 1.5 km s^^. A more accurate value has 
been determined from spectra collected with the UVES spectro- 
graph: wsini = 13.0 ± l.Okms"^ (Vladimir Elkin, private com- 
munication). These UVES spectra have better resolution and higher 
S/N than previous studies, thus we use this value of v sin i to deter- 
mine the inclination to be i = 37 ± 4°. 

The rotation of a Cir was first measured indirectly by Kurtz 
et al. (1994) from the amplitude modulation of the principal os- 
cillation mode at 2442 pHz. From the observed frequency separa- 
tion of the two symmetrical side-lobes, Kurtz et al. (1994) deter- 
mined the period Prot = 4.463 ± 0.010 d. This was refined to get 
Ptoi = 4.4790 ±0.0001 using photometry from several observing 
campaigns with a baseline of 12 y (Kurtz et al. 1994). These results 
are in very good agreement with our direct measurement. 

We interpret the modulation of the light curves in Fig. 2 
as being due to surface features, most likely two spots of over- 
abundances of rare earth elements located near the equator, thus 
not associated with the magnetic and pulsation poles. Such spots 
are typical of the CVn class of photometrically variable mag- 
netic Ap stars (see, e.g., Pyper 1969). Redistribution of flux from 
blue to red because of the increased line blanketing in the blue 
often gives rise to rotational light variations that are out-of-phase 
between blue- and red-filter observations in Ap stars. We cannot 
determine with certainty whether this is the case for a Cir, as we 
have only very broad-band observations with WIRE with an uncer- 
tain filter function (see section 3.4). Since the flux from the star is 
dominated by blue light, given that T^a = 7420 ± 170 K (Paper I), 
it is likely that the double wave variation shows minima when each 
spot is closest to the line-of-sight. The two spots are of different 
size, hence giving the double wave variation with rotation. We note 
that the double wave is asymmetrical (as can be seen in Fig. 2), 
since the secondary minimum lies at phase 0.519 ± 0.008. It is 
worth noting that the configuration of the spots that give rise to the 
rotational light modulation did not appear to change significantly 
during any of the runs with WIRE. We will discuss the spot config- 
uration in more detail in a future publication concerning Doppler 
imaging of the rotational light curve. 

Measurements of the magnetic field of a Cir are inconclusive 
regarding its variability with the rotation period. Bychkov et al. 



Table 3. Frequencies extracted from the two combined WIRE data sets from 
2006. Three frequencies were fitted to the SAAO _B-filter photometry as 
described in the text. The phases are with respect to to = HJD245 3939. 







WIRE 




SAAO 




f 


^WIRE 


'/'WIRE 


Ab 


<t>B 


ID 


[^Hz] 


[mmag] 


[0..1] 


[mmag] 


[0..1] 


h 


2442.0566(2) 


0.652(5) 


0.506(2) 


1.54(3) 


0.490(3) 


/2 


2265.43(2) 










h 


2341.79(1) 










U 


2366.47(4) 


0.047(8) 


0.67(4) 






h 


2566.98(3) 


0.044(7) 


0.18(3) 






h 


2411.8822(9) 


0.149(5) 


0.471(6) 


0.28(3) 


0.46(2) 


h 


2472.2272(5) 


0.188(5) 


0.041(4) 


0.42(3) 


0.06(2) 


fi 


2439.47(2) 


0.089(5) 


0.06(1) 






ft 


2444.64(2) 


0.079(5) 


0.97(2) 







Notes: /2 and /s were detected by Kurtz et al. (1994), but are not 
present in the WIRE data, /s was only detected in the Jul' 06 data. 

and correspond to /i ± /rot , where /rot is the rotational 
frequency. 



(2005) reanalysed sparse data collected from the literature span- 
ning 14 y and showed a rotation curve consistent with Prot = 
4.4794 d. We have repeated that exercise including a new measure- 
ment (Hubrig et al. 2004), making the time span 23 y, and find a 
marginal indication of magnetic variability with the known rota- 
tion period Prot = 4.4794 d, at the 2.5a significance level, if we 
neglect one outlying positive field measurement. This does not give 
us confidence that magnetic variability with the rotation period has 
actually been detected, and it is certainly not possible to give a time 
of magnetic maximum that can be compared to the time of rota- 
tional light minimum, or to the time of pulsation maximum. 

We can conclude that the longitudinal field for a Cir has been 
detected. The mean value derived from the data of Hubrig et al. 
(2004) is the most reliable at -239 ± 27G (internal error). The 
evidence from the data in the literature is that the longitudinal field 
is probably not reversing; i.e. only one magnetic pole is seen over 
the rotation period, implying that i + /3 < 90°, where i is the 
rotational inclination and 13 is the magnetic obliquity of the negative 
pole. For oblique magnetic rotators, the strength of the measured 
longitudinal magnetic field is proportional to cos a, where a in the 
angle between the magnetic pole and the line-of-sight. That varies 
as 

cos a = cos i cos /3 + sin i sin /3 cos fit, (1) 

where is the rotational frequency. 

Our rotational phase zero - the time of rotational light mini- 
mum in Fig. 2 - is therefore the time when the largest of the two 
equatorial spots is on the observed hemisphere closest to the line- 
of-sight. Since two equatorial spots are needed to describe the light 
curve, they cannot be associated with the one magnetic pole that is 
always visible. The pulsation pole is also not associated with the 
spots, as we will see in Sect. 4. This is not normal for an roAp star, 
hence for all interpretation related to the rotational light curve, we 
proceed with caution. 

3.2 Analysis of the rapid oscillations 

We analysed the rapid oscillations in the WIRE and SAAO light 
curves using the Fourier analysis package PERIOD04 (Lenz & 
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Figure 3. Prewhitening of the WIRE and SAAO data sets. In each row the triangles and veilical hnes mark the frequency being subtracted and the hatched grey 
regions mark the harmonics of the WIRE orbit. In the amplitude spectra from 2000 and 2005 we detect only the dominant /i mode. All four amplitude spectra 
from 2006 show /i, /g and fi, although for the SAAO data the S/N is poor for the two weaker modes. In the bottom row four additional peaks, /i ± frot, 
/4, and /s, are marked in the WIRE data sets from 2006. 



Breger 2005) and a Discrete Fourier Transform (Kurtz 1985), and 
we found excellent agreement between the two methods. We first 
subtracted the rotation signal, as found above. We initially anal- 
ysed each data set independently, and the prewhitening process is 
shown in Fig. 3. The top panels show the amplitude spectra of 
the six observed light curves, i.e. four from WIRE and two from 
SAAO. In each row of panels the frequency with the highest am- 
plitude is marked by the triangle and a vertical line. The hatched 
regions mark the harmonics of the orbital frequency of WIRE, 
which changed considerably from Sep'OO (/wire = 174.2 pHz) 
to Jur06 (/wire = 178.9 pHz) due to the decay of the orbit. The 
signal seen at the orbital harmonics (n x /wire) has leaked from 
drift noise at low frequencies. Note that in the Sep'OO spectrum, the 
main /i mode lies close to an orbital harmonic (?i = 14). The sec- 
ond row of panels in Fig. 3 shows the residual amplitude spectra 
after subtracting the /i mode at 2442 /iHz. No additional signifi- 
cant peaks are detected in the WIRE data from Sep'OO and Feb'05. 
However, both the WIRE data sets from 2006 and the simultane- 
ous SAAO data show two peaks, /e and fr, that have not been 
detected in previous campaigns. Interestingly, these peaks lie sym- 
metrically at ±30 /iHz around the main mode. Note that in the am- 
plitude spectra from SAAO, the 1 cd^^ alias peaks are sometimes 
slightly higher than the "true peak" (known from the WIRE data and 
marked by the vertical line). This problem is well known for single- 
site data, especially at low S/N. Fortunately, there is no ambiguity 
in our case since we have the high-S/N satellite data. In the bottom 
panels in Fig. 3 we see no signal above the noise in either of the 
SAAO data sets. In the Feb' 06 and Jul' 06 WIRE data sets we detect 



additional peaks^ at f^=fii: /rot and /4 = 2366 /iHz, while in 
Jur06 we detect also a low amplitude peak at /s — 2566 /tHz. 



3.3 Combining the 2006 WIRE data sets 

Since the two WIRE data sets from 2006 contain nearly the same 
peaks in the amplitude spectrum, we combined the data sets to be 
able to extract more precise frequencies: the formal error on fre- 
quency scales with l/Tobs, where Tobs is the baseline of the obser- 
vations (e.g. Montgomery & O'Donoghue 1999). There is a caveat 
when combining these data due to the 4-month gap between them, 
resulting in a spectral window with many sub-peaks that have al- 
most equal amplitude. To illustrate the problem, spectral windows 
are shown in Fig. 4 for the Feb'06 data set and the combined Feb'06 
+ Jul' 06 data sets. In the top panel the detected frequencies are 
marked by arrows, including /2 and /a, which were observed by 
Kurtz et al. (1994) but are not present in the WIRE data sets. 

To see if there is any ambiguity when extracting the frequen- 
cies, we repeated the analysis for 100 simulations of the combined 
Feb'06 + Jur06 time series, to see how often the input frequencies 
were correctly identified. To make the simulations, we followed 
the approach of Bruntt et al. (2007, their Appendix B). The simula- 
tions include the observed frequencies, the rotational double wave, 
white noise, and a 1// noise component consistent with the obser- 
vations. The simulations were analysed by an automatic prewhiten- 
ing procedure. The inserted frequencies were retrieved in all 100 
cases, except 11/100 for fi and 3/100 for (these modes have 
S /N ~ 6 and 10). In these cases frequencies offset by ±0.083 pUz 



1 /rot = 1/Prot = 2.5841 ± 0.0001 AtHz. 
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Frequency [/xHz] 

Figure 4. Spectral windows for the Feb'06 data set {top panel) and the combined Feb'06 + Jul'06 data sets {bottom panel) from WIRE. The insets show details 
of the central peak. The arrows mark the positions of the frequencies detected in a Cir and the hatched regions mark the location of the orbital harmonics. 
Note that /2 (close to an alias peak of /i) and /a are from Kurtz et al. (1994) and were not detected in the WIRE data sets. 



were picked by the automatic procedure, corresponding to an alias 
peak offset by il/Tgap, where Tgap is the duration of the gap. Due 
to these relatively high false-alarm ratios, we will quote results for 
the low amplitude modes using only the best data set from Jur06. 

In Table 3 we give the final results for the combined data set 
for /i, /e and /t, while the results for the other modes are based on 
the Jur06 data set. Note that the frequencies were fixed, with /i 
from the combined data set and /rot = l/^rot from Sect. 3.1. We 
used the results of the simulations to determine the uncertainties of 
the frequency, amplitude, and phase calculated as the standard de- 
viation of the extracted parameters. The uncertainties are given in 
the parentheses in Table 3, and are in general agreement with the- 
oretical expressions, assuming pure white noise (e.g. Montgomery 
& O'Donoghue 1999). 



/e and /y are in agreement with this shift, but the uncertainties are 
three times larger and cannot be used to constrain the phase shift. 

Kurtz & Balona (1984) made simultaneous observations in 
Johnson B and V and measured only small phase shifts <j>v — 
(pB-v = +13° ±6° and (j!>v — 0s = — 7M± 5? 1. They measured 
the amplitude ratio for /i to be Ab/Av = 2.28 ± 0.26, which is 
essentially the same as our ratio of Ab /Awire- Thus the ampli- 
tude in the WIRE data is similar to that which would be observed 
through a Johnson V filter, strongly suggesting that the bandpass 
of the WIRE photometry is centered at about the same wavelength 
as Johnson V. The interpretation of the small phase shifts is more 
problematic, since it is well-known from high resolution spectro- 
scopic studies that the pulsation phase in roAp stars can be a sen- 
sitive function of atmospheric depth of the observations (see, e.g., 
Kurtz et al. 2006a, Ryabchikova et al. 2007). 



3.4 Calibration of the WIRE photometry 

To calibrate the measured amplitude and phases of the WIRE pho- 
tometry, we analysed the combined ground-based Johnson B data 
from the two SAAO runs, obtained at the same time as the two 
WIRE data sets from 2006. We made simulations of the data as de- 
scribed in Sect 3.3, which showed that the automatic prewhitening 
was seriously affected by the complicated spectral window with 
strong Icd^^ alias peaks. Fortunately, the frequencies of /i, fa 
and fr are determined very accurately from the WIRE data and we 
therefore held the frequencies fixed when determining their ampli- 
tudes and phases. The fitted values are listed in the last two columns 
in Table 3. 

We measured the amplitude ratio and the phase shift between 
the WIRE and SAAO observations. The amplitude ratios for the 
three modes /i, fe and /y are: Ab /Awire = 2.36 ± 0.11, 1.88 ± 
0.38, and 2.23 ± 0.36, with a weighted mean {Ab /Awire) = 
2.32 ± 0.10. The results for the /i mode shows a small but signif- 
icant phase shift of ^wire — <^b ~ +5?8 ± 1?3. The phases of 



4 IDENTIFICATION OF THE PRIMARY MODE 

In the following we will discuss the possible mode identification of 
the primary /i mode using the relative amplitudes of the rotational 
sidelobes. This is done in the framework of the oblique pulsator 
model, which postulates that the axis of the non-radial pulsation 
modes is inclined to the rotation axis of the star. In some versions 
of the model the pulsation axis coincides with the magnetic axis 
(Kurtz 1982, Shibahashi 1986, Kurtz & Shibahashi 1986, Shiba- 
hashi & Takata 1993, Takata & Shibahashi 1995), while in another 
version the pulsational axis may be inclined to both the rotational 
and magnetic axes (Bigot & Dziembowski 2002). Here, we inter- 
pret the rotational frequency triplet /i ± /rot in terms of the for- 
malism of Shibahashi (1986) and Kurtz & Shibahashi (1986). 

We expect the observed modes in photometry to be of low 
degree (typically I ^ 2), given that generally perturbations associ- 
ated with modes of higher degree cancel out when averaged over 
the stellar disk. However, this may not always be the case in roAp 
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Figure 5. The top panel is the Hght curve of a Cir from WIRE in Jul' 06. 
The middle and bottom panel shows variation of the amplitude and phase 
of the primary mode at 2442 /iHz determined by fits to subsets of the light 
curve. The vertical lines mark the time zero points at maximum amplitude 
(tmax; dashed line) and minimum light (trot ; solid line). 



Stars, since the presence of strong magnetic fields that permeate the 
surface of these stars may resuh in a distortion of the eigenfunc- 
tions at the surface (Dziembowski & Goode 1996; Cunha & Gough 
2000; Saio & Gautschy 2004). Thus, a mode of degree higher than 
I — 2 could in principle give rise to lower-degree components at the 
surface and still be detected. If the magnetic perturbations are sym- 
metric about the magnetic equator, as is the case for pure dipolar or 
pure quadrupolar magnetic topologies, the additional components 
of the perturbed eigenfunctions will have the same parity (odd or 
even as the non-perturbed one. We will discuss this further in the 
third paper of this series (Brandao et al., in preparation), in which 
detailed magnetic modelling of the pulsations in a Cir will also be 
presented. 



4.1 Properties of the rotational sidelobes 

We observed two frequency sidelobes around the principal pul- 
sation frequency /i, separated by 2.58 /iHz, which is the same 
as the rotational frequency. In the two amplitude spectra from 
WIRE 2006 in Fig. 3 the sidelobes at /i ± /rot are clearly seen. 
To measure the amplitudes and phases of the sidelobes we as- 
sume they are separated by exactly the rotation frequency, as re- 
quired by the oblique pulsator model. We use the highest pre- 
cision value we have for the rotational frequency (cf. Sect. 3.1; 
/rot = 2.5841 ± 0.0001 fiHz) and fitted an equally spaced fre- 
quency triplet about /i = 2442.0566 /^Hz by least squares, with 
the requirement that the zero point of the time scale be chosen such 
that the phases of the rotational sidelobes be equal (i.e. (j)^ — (j>i), 
as expected for axisymmetric modes in the oblique pulsator model 
(Kurtz 1982). 

The results of that fit are shown in Table 4 for the zero point 
in time tmax = HJD245 3935.91562. The amplitudes of the 
two rotational sidelobes are equal to within 1.5 a, i.e. Af /A^ = 
0.854±0.096. This is consistent with the magnetic perturbations to 
the eigenfrequencies being very much stronger than the rotational 
perturbations. The zero point in time for the least-squares fit was 
chosen such that the phases of the rotational sidelobes are equal, 
but the phase of the central frequency of the triplet was not con- 
strained by this choice of fmax. The nearly identical phases (they 
differ by 1.6 a) are thus consistent with a purely geometrical am- 



Table 4. A least squares fit to all of the WIRE 2006 data of the principal 
frequency, fi = 2442.0566 fiHz and its two rotational sidelobes which - 
within the oblique pulsator model - are separated by exactly the rotation 
frequency of /rot = 2.5841 ± 0.0001 ^tHz. 
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A 


4> 


ID 


fiHz 


mmag 


[0..1] 


/r 


2439.4725 


0.082 ± 0.006 


0.423 ±0.011 


/i 


2442.0566 


0.618 ±0.006 


0.4050±0.0014 


ft 


2444.6407 


0.070 ± 0.006 


0.423 ±0.013 



plitude modulation with rotational aspect in the oblique pulsator 
model. 

To investigate the variation of the amplitude we followed the 
approach of Kurtz et al. (1994) to analyse the WIRE Jur06 data 
set, which has the highest S/N of all of our data. We first sub- 
tracted the other modes from the data set, i.e. the rotational dou- 
ble wave /a + /s (Sect. 3.1) and the low amplitude modes fA-fr 
(Sect. 3.2). We then split the data into subsets containing 30 pulsa- 
tion cycles of /i , corresponding to time bins of Af ~ 3.4 h. In each 
subset we found the amplitude and phase of the /i mode by least 
squares, holding the frequency fixed. We find that the amplitude of 
/i is modulated by the rotation while the phase is not. The result 
is shown in Fig. 5, where the variations in amplitude and phase of 
/i are shown in the middle and bottom panels, respectively. The 
top panel shows the WIRE light curve for comparison. The time 
of minimum light (trot, as defined in Sect. 3.1) is marked with a 
solid line and the time of maximum amplitude (tmax) is marked 
with a dashed line. The delay in the maximum amplitude with re- 
spect to the zero point for the rotation phase (Sect. 3.1) is = 
^max — 0rot = 0.287±0.003. In Other words, the time of maximum 
amplitude occurs at rotation phase 1 — A(f> = 0.713 (see Fig. 2). 
Comparison of the time zero points for the rotational light varia- 
tion and the pulsation maximum gives nearly the same result: the 
rotational phase is ^irot = (imax — irot)/frot ~ —0.289 = 0.711. 

For an axisymmetric dipole mode (/ = 1, m = 0) aligned with 
the magnetic field, and for spots centred on the magnetic poles, we 
expect the time of pulsation maximum to coincide with the time of 
rotational light minimum, whereas, as can be seen in Fig. 5, rota- 
tional phase 0.71 coincides with one of the rotational light maxima. 
The essentially pure amplitude modulation with a sinusoidal varia- 
tion over the rotation cycle is consistent with a pure dipole pulsation 
mode, leading to the conclusion that the spots are not concentric 
with the magnetic field. While small departures between the spots 
and magnetic poles do occur in some magnetic Ap stars, e.g. in the 
roAp star HR3831 (Kurtz et al. 1992), the case of a Cir requires 
equatorial spots far from the pulsation pole. For further progress, 
new, very precise measurements of the magnetic field of a Cir are 
needed to show a convincing rotational variation, as we pointed out 
in Sect. 3.1. 



4.2 Is /i a dipole or a quadrupole mode? 

The observed pulsation amplitude of an obliquely pulsating nonra- 
dial mode varies with rotation, as the aspect of the mode changes. 
This gives rise to a frequency multiplet, separated by the rotation 
frequency, that describes the amplitude modulation (and also phase 
modulation, if present). Within the formalism of the oblique pul- 
sator model (Shibahashi 1986; Kurtz & Shibahashi 1986; Kurtz 
et al. 1990) we have calculated the rotational amplitude modula- 
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Figure 6. The relative amplitude modulation over the rotation cycle for two 
dipole modes, / = 1, m = (solid circles) for /3 = 18° and I = l,m = 
1 (open circles) for /3 = 108°, and the axisymmetric quadrupole mode 
with i = 2,m = (crosses) for f3 = 51° (in all cases i = 37°). The 
values of /3 have been derived from the constraints of the amplitude ratio 
of the rotational sidelobes to the amplitude of the principal frequency as 
described in detail in the text. For each mode the maximum amplitude has 
been normalised to 1 . 



tion expected for three possible pulsation modes, = 1, m = 0), 
{I = l,m = ±1), and {I = 2,m = 0), such that the amplitude 
modulation agrees with the observed amplitudes of the rotationally 
split frequency triplet given in Table 4. The results are shown in 
Fig. 6 and each case will be discussed in detail below. 



4.2.1 An axisymmetric dipole mode (I = l,m = Q) 

Whatever the relationship between the rotational light curve and 
the pulsation, the frequency triplet of Table 4 is consistent with an 
axisymmetric dipole mode, or any mode of higher degree to which 
distortion by the magnetic field could result in the addition of a 
significant {I = l,m = Q) component. This leads to the constraint 
(Kurtz et al. 1990): 



— = tan i tan = 0.246 ± 0.014 = n , 



(2) 



where i is the rotational inclination and /3 is the magnetic obliquity 
of the negative pole, and A+i, A-\ and refer to the amplitudes 
of the rotational sidelobes and the principal frequency. Since we 
have already determined the rotational inclination to be i = 37±4° 
in Sect. 3.1, that leads to /3 = 18 ± 3°. These values are consistent 
with only one magnetic pole being visible over the rotation cycle, 
i.e. with i + P < 90° as we have discussed in Sect. 3.1. Fig. 6 
shows an axisymmetric dipole mode with these values of /3 and i 
(solid circles). It has a sinusoidal amplitude modulation, which is 
in agreement with the observed modulation in Fig. 5. 

The oblique rotator model for a dipole magnetic field leads to 
a similar constraint on the geometry; Hubrig et al. (2007) found 
tan i tan/? = 0.66 with no estimated error. We cannot, therefore, 
assess whether there is an inconsistency between these derivations 
of tan i tan /9, but note that Hubrig et al. (2007) find i = 38°, 
/3 — 40°, hence i+/3 < 90°, also indicating that only one magnetic 
pole is visible. 



A+i+A-i 
Ao 



tan I 
tan P 



= n, 



(3) 



which leads to /? = 108 ± 3° for i = 37 ± 4°. This implies that 
both magnetic poles should be observed, contrary to observations. 
The amplitude modulation of a sectoral I = l,m — 1 mode with 
these values of /3 and i is shown in Fig. 6 (open circles) and exhibit 
a double wave that is not observed; this is because both poles are 
seen, but from different aspect. The pulsation mode can therefore 
not be a dipole mode with m = ±1. 

4.2.3 An axisymmetric quadrupole mode I = 2,m — 

From our theoretical modelling in Sect. 5.1 we show that both even 
and odd I modes are present, so it is worth exploring whether the 
principal mode could have even degree, seen by the observer as 
a quadrupole mode. For an axisymmetric quadrupole mode with 
Z = 2, m = we expect a frequency quintuplet with separations of 
/rot ~ 2.5841 fjHz. Since the outer two frequencies of this quin- 
tuplet are not observed, we assume that their amplitudes are below 
the 3ct detection limit (i.e. A±2 < 0.018 ± 0.002 mmag) and solve 
for the magnetic obliquity that is consistent with this. We then ex- 
pect (Kurtz et al. 1990; Shibahashi 1986) 

A-f-i + A-i 12 sin /3 cos /3 sin i cos i 



Ao 



(3cos2/3- l)(3c 



1) 



ri, 



and 



r2 



A+2+A. 
Ao 



3 sin^ (3 sin^ i 



(3cos2/?-l)(3cos2i-l) 
s; 0.058 ± 0.005. 

We can divide Eq. (6) by Eq. (4) to yield 



(4) 

(5) 

(6) 
(7) 

(8) 



4r9 

tan i tan/? = — ^ 0.94 ±0.10. 

ri 

For i = 37 ± 4°, Eq. (8) yields /3 = 51 ± 5°. This also fulfills 
the requirement that i + l3 < 90°, but the amplitude variation is 
clearly not sinusoidal with rotation, as can be seen in Fig. 6 (cross 
symbols). In fact, for this quadrupole the amplitude goes through 
zero and reverses phase over the rotation cycle, in clear conflict 
with the observations. We therefore conclude that the principal fre- 
quency does not have the observational properties of an axisym- 
metric quadrupole mode. 



4.3 Concluding remarks on the mode identification of /i 

The apparent lack of any rotational phase variation (Sect. 4.1) and 
the rotational amplitude variation of the principal pulsation mode 
(Sect. 4.2) strongly support the supposition that the principal mode 
has the observational properties of an axisymmetric I = l,m = 
dipole mode. Note well, however, that we do not constrain our theo- 
retical modelling in the next section to dipole modes for the princi- 
pal pulsation mode, since magnetic perturbations may make modes 
of a different degree appear to be very similar to dipole modes ob- 
servationally. 



4.2.2 A non-axisymmetric dipole mode (I = l,m = ±1) 

Non-axisymmetric dipole modes (I — l,m = ±1) within the 
oblique pulsator model constraints on the geometry (Kurtz & 
Shibahashi 1986, Shibahashi 1986) give 



5 ASTEROSEISMOLOGY OF a Cir 

In the following we will compare theoretical pulsation models with 
the new observational constraints consisting of the improved fun- 
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Figure 7. Frequencies detected in a Cir from tlie WIRE data and /2 / /a 
from Kurtz et al. (1994). The vertical dashed lines mark the mean separation 
between /g, /i and /y. 



damental parameters (Paper I) and the detection of the large sepa- 
ration in the current work. 

Before comparing the observed and theoretical frequencies 
it is crucial to include only the observed frequencies for which 
there are no ambiguities about their detection. In Fig. 7 we show 
a schematic frequency spectrum of seven modes that we believe 
are due to pulsation and intrinsic to a Cir. In our list of secure fre- 
quencies in Table 3 we have included the five frequencies /1-/5 
detected by Kurtz et al. (1994) and the two new frequencies /e and 
fr. In Fig. 7 we have scaled the amplitudes to the B filter system 
using the ratio Ab /Awire from Sect. 3.4, and for /2 and /a we 
use the scaling Av/Ab = 1.29 (Kurtz et al. 1994). The vertical 
dashed lines correspond to the average separation between consec- 
utive modes in the triplet fa + fi + /y. 

The frequencies (2 and /a were not detected in our data al- 
though from the observed noise level in the WIRE 2006 data and 
the observed amplitude we should have detected them with a S/N 
ratio around 6-8. A likely explanation is that their amplitudes have 
changed over time. Alternatively, it is possible that /2 was not de- 
tected due to its proximity to the sideband f2 — fi — /wire 
(A/ = 2.3 /iHz), while the detection of /a was perhaps not possi- 
ble due to it being relatively close to an orbital harmonic (n — 13, 
A/ = 16.2 /iHz). These possibilities are illustrated in Fig. 4. 

The frequencies /e and fr must have appeared some time 
between Feb'05 and Feb'06, since they were not detected in the 
WIRE dala sets from Sep'OO and Feb'05. If their amplitudes had not 
changed, the peaks should have been detected, with S/N of about 
14 and 7 in the amplitude spectra from Sep'OO and Feb'05, respec- 
tively. Nevertheless, as they are seen in both of the simultaneous 
runs with WIRE and from SAAO in 2006, they must be intrinsic to 
Q Cir. 

A few additional frequencies have been reported based on 
spectroscopic campaigns (Balona & Laney 2003; Kurtz et al. 
2006b), but we have not included them in the modelling. The fre- 
quencies reported by Balona & Laney (2003) are clearly affected 
by 1 cd~^ aliasing problems, and the data analysis by Kurtz et al. 
(2006b) is complicated by the relatively poor frequency resolution 
while at the same time the frequencies lie close to the principal 
mode. 

The frequencies shown in Fig. 7 contain information about the 
average properties of a Cir as well as details about its internal struc- 
ture and magnetic field. In the following we will attempt a forward 
comparison with the theoretical oscillation spectra obtained from 
models. To avoid confusion with the observed frequencies /, we 



shall designate the theoretical frequencies by Vnj, where n and I 
are, respectively, the radial order and the degree of the eigenmodes. 



5.1 The large separation 

The most commonly used observable quantity when studying high- 
overtone p mode oscillations in stars is the large frequency separa- 
tion Ai/n.i, defined as the difference between modes of the same 
degree and consecutive radial orders: Ai^„,i = — i^n,i- The 

large separation is a function of frequency, but for high radial or- 
ders the dependence is small. In fact, in a spherically symmetric 
star, the frequencies of linear adiabatic acoustic pulsations of high 
radial order are expected to follow the asymptotic relation (Tassoul 
1980): 

~ Ai'o(n + V2 + l/4 + a) + en,i, (9) 

where Az^o ~ (2 J^^ dr/c)^^, i? is the stellar radius, c is the sound 
speed and r is the radial coordinate. The value of a depends on the 
properties of the reflection layer near the surface, and is a small 
term that depends mostly on the conditions in the stellar core. It can 
be seen that Ai/q is the inverse of the sound travel time across the 
diameter of the star. 

The high frequency of the oscillations observed in a Cir in- 
dicates that these are indeed acoustic oscillations with high radial 
order. Thus one may expect to find in the oscillation spectrum a 
regular pattern corresponding to the asymptotic relation. However, 
the asymptotic analysis assumes spherical symmetry of the star In 
a Cir, as in other roAp stars, the presence of a strong magnetic 
field introduces a source of deviation from spherical symmetry that 
can shift the oscillation frequencies significantly from the asymp- 
totic values (Cunha & Gough 2000; Cunha 2006; Saio & Gautschy 
2004). According to these studies, one may expect to still see regu- 
lar frequency spacings in some parts of the oscillation spectra, but 
abnormal spacing might also be found as in the well-studied case 
of HR 1217 (Kurtz et al. 2005; Cunha 2001). 

The observed frequencies fa and /y lie almost equidistantly 
around /i, to within a few nHz: /i - /e = 30.1746 ± 0.0009 /iHz 
and /t - ,/i = 30.1707 ± 0.0005 /iHz. We determine the aver- 
age separation from the slope of a first order polynomial fitted to 
the three frequencies: A/obs ~ 30.173 ± 0.004 /tHz. To be con- 
servative we adopt the difference between the two values as the 
uncertainty. This separation is marked by the vertical dashed lines 
in Fig. 7. We note that the frequencies /2 and /s appear to roughly 
agree with the spacing. 

If the regular frequency spacing of /e + /i + /t involves 
modes with degrees of different parity, and if their frequencies are 
not significantly modified by the magnetic field, then we may con- 
clude that the large separation is Avo ~ 2 A/obs- Otherwise, if 
only modes with degrees of the same parity are present, we have 
Ai^o ~ A/obs- Since Avo depends essentially on the mean density 
of the star, these two possibilities correspond to stars with very dif- 
ferent mean densities. It is straightforward to determine, based on 
non-asteroseismic observables of a Cir, which one of these possi- 
bilities is more plausible. 



5.2 The pulsation model 

To compute theoretical pulsation frequencies we first calculated a 
reference equilibrium model of a Cir using the stellar evolution 
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Table 5. Global parameters of the CESAM model used to represent a Cir. 
As input parameters we used = 0.70, Vq = 0.28, o = 1.6 and no 
overshooting. and Iq are, respectively, the initial hydrogen and helium 
abundances and a is the mixing length parameter The mass, M, luminos- 
ity, L, and radius, R, are expressed in solar units. is the homologous 
constant defined in Eq. (10). 



Af/M0 


L/Lq 


Teff [K] RIRq 


[/iHz] 


1.715 


1.019 


7421 1.96 


130.8 



code CESAM^ (Morel 1997). The fundamental parameters used for 
this model are given in Table 5 and are adopted from our analysis 
of interferometric data and spectra calibrated in flux as described 
in Paper I. The position corresponding to these parameters is given 
by the asterisk in the Hertzsprung-Russell (HR) diagram in Fig. 8. 

To calculate the theoretical oscillation frequencies for this 
model we used the linear adiabatic oscillation code ADIPLS 
(Christensen-Dalsgaard 2008). Since the observed frequencies are 
above the acoustic cutoff frequency for the reference model, we 
have applied at the outermost point of the pulsation models a fully 
reflective boundary condition of the type 5p — 0, where Sp is the 
Lagrangian perturbation to the pressure. The physical mechanism 
responsible for the reflection of the oscillations in roAp stars has 
been under debate for decades. It is likely that the reflection of 
the oscillations results from the direct effect of the magnetic field 
(Sousa & Cunha 2008). Despite the small value of the longitudi- 
nal magnetic field of a Cir (Hubrig et al. 2004; cf. Sec. 3.1), we 
note that the magnetic field intensity of this star is not necessarily 
small. If we assume a dipolar topology and the magnetic obliquity 
and rotational incUnation derived in Sect. 4.2.1, we find a polar 
magnetic strength of ~ 1 kG. Moreover, if the magnetic field has a 
quadrupolar component, its polar intensity will be even stronger A 
comparison between magnetic and gas pressures in the outer lay- 
ers of a magnetic stellar model with parameters similar to a Cir 
was presented by Cunha (2007). From their discussion it is clear 
that even a 1 kG magnetic field has a strong effect on the dynam- 
ics of the oscillations throughout the whole atmosphere of a typical 
roAp star Thus, it is expected that through mode conversion part 
of the mode energy is passed onto magnetic waves in the atmo- 
sphere. The energy transferred to these magnetic waves is kept in 
the mode. Thus, if the energy input in each cycle through the exci- 
tation mechanism (Balmforth et al. 2001; Cunha 2002; Saio 2005) 
is enough to compensate for the fraction of the mode energy that 
is lost through running acoustic waves in the atmosphere (Sousa & 
Cunha 2008) and through running magnetic waves in the interior 
(Roberts & Soward 1983), the mode will be over-stable and may 
reach a detectable amplitude. 

The use of a fully reflective boundary condition is an artifi- 
cial way to induce reflection of pulsations in the theoretical model 
and may shift the computed frequencies of individual modes away 
from their true values (Cunha 2006). Nevertheless, its effect on 
the individual mode frequencies is expected to be systematic and, 
thus, should cancel out when computing pulsating quantities that 
are based on frequency differences, particularly when these differ- 
ences involve only modes of the same degree. Since this is indeed 
the case for the large separation, our computed values will not be 
significantly affected by the adopted boundary condition. 

2 The CESAM code is available at 

http://www.astro.up.pt/corot/models/cesam/ 
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Figure 8. H-R diagram with three evolution tracks from CESAM. The po- 
sition of our reference model is marked with an asterisk and four additional 
models considered in our analysis are labeled A, B, C and D. The rectangle 
marks the 1 u uncertainty on T^ff and L/\jq and the two diagonal lines 
correspond to constant radii consistent with the 1 a uncertainty on the in- 
terferometric measurements from Paper I. 



The average large separation of the reference model is 
(Ai^„,i) = 62.5 /iHz. Modes of degrees / = 0, 1, 2 and 3, with 
frequencies around the observed values, were used to calculate 
this average. This value is close to 2A/obs and so the solution 
Avo ~ A/obs is firmly rejected. This means that modes with de- 
grees of alternating parity must be involved in the regular frequency 
pattern observed in the triplet /e, /i, /?■ 

It is beyond the scope of the present work to make a de- 
tailed model of a Cir. However, it is worth inspecting how the re- 
sults for our model compare with those for similar models placed 
within the region of uncertainty associated with the global pa- 
rameters of a Cir. We have considered four additional models, 
labeled A, B, C and D in Fig. 8. These models lie at the ex- 
tremes of the 1 cr-uncertainty rectangle (T^s — 7420 ± 170 K, 
L/Lq = 10.51 ± 0.60) and inside the strip defined by the two 
lines of constant radius corresponding to R/Kq = 1.967 ± 0.066 
(the interferometric radius from Paper I). To estimate the large sep- 
arations for these four models, we used the fact that the averaged 
asymptotic large separations follow a homologous relation of the 
type 

{Au„,i) = a ■ {R/Rq)-^^\ (10) 

The value of C',^ for our reference model of a Cir is 130.8 /xHz. In 
Table 6 we list the masses of these models, as derived from CESAM 
evolution tracks using the same input parameters as for our refer- 
ence model, their radii, as derived from their luminosities and effec- 
tive temperatures, and the corresponding large separations, derived 
through the homologous scaling in Eq. (10). 

Inspection of Table 6 shows that the large separations of 
the models (59 to 65 /xHz) agree well with the observed value 
(60 /iHz). However, this is not the case for the high effective tem- 
peratures that have been suggested in the literature (see Paper I for 
a discussion). For an effective temperature of 7900 K, and with 
a similar luminosity and mass, using the homologous scaling in 
Eq. (10) we find that the large separation of the model is about 
75/iHz, greatly exceeding the observed value. Hence, the astero- 
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Table 6. Masses, radii and average large separations obtained from the ho- 
mologous scaling given for models A, B, C, and D shown in Fig. 8. 
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Figure 9. The frequency differences defined by Eq. (11) are shown ver- 
sus mode frequency. Results ai'e given for the cases when the triplet is as- 
sumed to involve modes of degrees / = 1, 0, 1 (plus symbols), / = 2, 1, 2 
(crosses), i = 3, 2, 3 (asterisks), and / = 3, 0, 3 (diamonds). The observed 
value is marked by the horizontal dashed line and the hatched region marks 
the location of the triplet. 



seismic modelling of the observed large separation strongly sup- 
ports the effective temperature derived for a Cir in Paper I. 



5.3 Evidence for magnetic perturbation of the modes 

Despite the general agreement found between the observed and 
computed large separations, there are aspects about the observed 
frequencies that cannot be accounted for with the models consid- 
ered here. Firstly, only the frequencies /e, /i and /y seem to follow 
a regular frequency pattern. In particular, the low amplitude mode 
fi, which is seen in both the WIRE data presented in this paper and 
the data acquired by Kurtz et al. (1994), is offset by about a quarter 
of the large separation from where it would be expected, according 
to the asymptotic approximation. Moreover, the almost-equal value 
of the separations (/i — /e) and (/y — /i) is also not reproduced 
by the reference model. If the triplet is assumed to correspond to 
modes of alternating parity and adjacent orders, then the difference 
between these two frequency separations can be written 



(11) 



Under the same conditions as those required for Eq. (9) to be valid, 
5vn,i is asymptotically given by (e.g. Cunha et al. 2007) 

Ai^o Ac dr 



5v, 



-4(/ + l). 



(12) 



Given the increase of the sound speed with depth in the star, Sv^.i 
is expected to be positive, increasing proportionally to (/ + 1). 

The observed value for this quantity is Suohs = ~ ((/? ^ fi)^ 
(/i - /e)) = 0-0039 ± 0.0014 /iHz. In Fig. 9 we plot Siy„,i for our 
reference model (cf. Table 5), assuming different combinations of 
the degrees of the modes in the triplet: I — (1, 0, 1), I — (2, 1, 2), 
I = (3, 2, 3), and I = (3, 0, 3). In the last case SPn,i is defined in a 



similar way, but with {n — 1,1 + 1) and {n, I + 1) replaced by (n — 
2,1 + 3) and (n— 1, Z + 3), respectively. Except for the combination 
/ — (1,0, 1), the values of Si/nj are positive and increase with 
as expected from the asymptotic expression. The deviation of 
5vn,i from the expected trend in the case of alternating modes of 
degrees / = 1, Z = and / = 1 is not surprising, since in the 
asymptotic analysis leading to Eqs. (9) and (12) it is assumed that 
the sound speed varies smoothly with depth. Because A-type stars 
have convective cores, that assumption breaks down at the edge of 
the core. The modes of lowest degree, which propagate deeper, will 
feel the effect of that rapid variation in the sound speed, and thus, 
their relative frequency is likely to deviate more strongly from the 
asymptotic behavior. 

Regardless of how closely the quantity Si'n,i derived from our 
model follows the asymptotic behavior expressed by Eq. (12), it is 
clear that for all possible combinations of modes considered, this 
quantity is significantly different from the observed value (dashed 
horizontal line in Fig. 9) in the region where the modes are observed 
(the hatched region from 2.4 to 2.5 mHz). The smallest absolute 
value found for the theoretical Siyn,i in the range of the observed 
frequencies, is around 2.0 /iHz for i = (1,0, 1), which is much 
larger than the 3.9 nHz determined from the observations. Consid- 
ering modes of degree higher than / = 3 will not solve this discrep- 
ancy, since Svn.i will increase with increasing I. Moreover, we find 
that this discrepancy is also present for models A, B, C, and D, i.e. 
none of our theoretical calculations can explain the almost-equal 
separation observed between the three frequencies fa, fi and /y. It 
remains to be investigated whether the effect of the magnetic field, 
which has been ignored in the models considered so far, may ex- 
plain this and the other apparent disagreements between theoretical 
and observed quantities discussed here. 

The effect of the magnetic field on the pulsations of roAp stars 
has been studied by different authors for over a decade (Dziem- 
bowski & Goode 1996; Bigot et al. 2000; Cunha & Gough 2000; 
Saio & Gautschy 2004; Cunha 2006; Sousa & Cunha 2008). How- 
ever, it is only recently that magnetic models have started to be 
used, in a forward approach, in an attempt to fit the asteroseismic 
data of particular pulsators (Gruberbauer et al. 2008; Huber et al. 
2008). The study of pulsations in the presence of a magnetic field 
is a complex one. In order to make it treatable, a number of as- 
sumptions and approximations have to be introduced. Thus, quan- 
titative comparisons between observed frequencies and frequencies 
derived from magnetic models require some caution. Nevertheless, 
there are a number of robust results derived from the magnetic mod- 
els that one may confidently use. Cunha (2005) argued that the ef- 
fect of a dipolar magnetic field is always to increase the large sep- 
aration by a small amount, that can vary from nearly zero to a few 
fjHz. Since our non-magnetic models reproduce well the observed 
large separation, we expect that in the case of a Cir the change in 
the large separation when including the magnetic field will be rel- 
atively small. On the other hand, a detailed magnetic modelling of 
a Cir is necessary to test whether the effects produced by a mag- 
netic field with properties consistent with the observational con- 
straints may remove the discrepancies found between our mod- 
els and the observed value of 5i'n,i, as well as explain the devia- 
tion from the asymptotic trend of some of the observed oscillation 
frequencies. Such analysis will be presented in a separate paper 
(Brandao et al., in preparation). 

Finally, we would like to emphasize that the almost-equal sep- 
aration of the modes that compose the main triplet raises a question 
that goes beyond the simple issue of whether a particular magnetic 
field may bring the model predictions into agreement with the ob- 
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servations. In fact, even if such a magnetic field is found, its effect 
will depend on the frequency and so a precise cancellation of the 
difference (/t — /i) — (/i — /e) may occur only in a very particu- 
lar range of frequencies. Therefore, the question of whether there is 
an underlying reason for the principal modes to be found precisely 
in that range of frequencies, must also be considered. In particular, 
model simulations that provide a measurement of how often this 
precise cancellation is expected to occur need to be carried out, and 
compared with the empirical evidence from observations of other 
multi-periodic roAp stars. 

5.4 The Hipparcos and asteroseismic parallaxes 

Prior to the WIRE observations presented in this paper, the pat- 
tern of observed frequencies observed in a Cir by Kurtz et al. 
(1994) suggested a large separation of 50 /iHz. Matthews et al. 
(1999) used this value, together with an effective temperature of 
Toft = 8000 K and a mass of M = 2.0 ± 0.5 Mq to calculate the 
expected "asteroseismic" parallax of a Cir. The value they deter- 
mined, TTast = 51.3 ± 1.2, was significantly smaller (Jo) than the 
parallax measurement by ///pparcoi- (7rHip-97 = 60.97±0.58 mas; 
Ferryman & ESA 1997). Moreover, a similar calculation for an- 
other 1 1 roAp stars indicated that the Hipparcos parallaxes were 
systematically larger than the parallaxes derived from the astero- 
seismic data of these roAp stars. The new, higher value for the large 
separation found in this paper, together with the lower mass and 
Tefl from Paper I, resolve the discrepancy found for a Cir. Using 
Eq. (4) from Matthews et al. (1999) we derive a luminosity of a Cir 
of L = 10.6 ± 1.5 Lq. We adopt a mass of M = 1.7 ± 0.2 Mq 
and Tcff = 7420 ± 170 K (both values from Paper I), a bolomet- 
ric correction of BC = 0.19 ± 0.08 (Brandao et al. 2007; av- 
erage of their two estimates), Mboi.Q = 4.75, V — 3.19, and 
E{B - \/)max = 0.06 ± 0.02 (from Matthews et al. 1999) to get 
7i"ast = 62.8 ± 5.3 mas. This is in good agreement with the par- 
allax from Hipparcos, which was recently updated to tthip-ot = 
60.36 ± 0.14 mas (van Leeuwen 2007). This resolves the problem 
of the low asteroseismic parallax for a Cir found by Matthews et al. 
(1999). Moreover, the lower effective temperature of a Cir found 
in Paper I conoborates one of the suggestions of Matthews et al. 
(1999), namely that the systematic discrepancy between the aster- 
oseismic and Hipparcos parallaxes could be associated with a sys- 
tematic error in the determination of the effective temperatures of 
roAp stars. 

In Fig. 10 we plot the Hipparcos parallaxes versus the aster- 
oseismic parallaxes for eight roAp stars. This is an update of the 
original Fig. 1 from Matthews et al. (1999). Four stars are omit- 
ted (HD 119027, 166473, 190290, and 203932) because the un- 
certainty on their Hipparcos parallax is 8-10 times larger than for 
the other stars. For a Cir we used the values found above, and 
for the other stars we used the information in Tables 1 and 2 of 
Matthews et al. (1999) and updated parallaxes from van Leeuwen 
(2007). Note that we have adopted larger uncertainties on T^ff 
(or = 200 K) and on mass (aM ~ 0.5 M©). The solid line in 
Fig. 10 shows equality, and it seems that there is now acceptable 
agreement for most stars, the exceptions being HD 217522 and 
HD 137949. 



6 CONCLUSIONS 

We have analysed observations of the roAp star a Cir, compris- 
ing 84 days of photometric time series from the star tracker on the 
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Figure 10. The Hipparcos pai'allaxes of eight roAp stai's are plotted vs. the 
asteroseismic parallaxes as determined from the large separation, T^b, and 
mass. This is an updated version of Fig. 1 from Matthews et af (1999). 



WIRE satellite. We used data from four different runs lasting be- 
tween 8 and 42 days collected over six years, or nearly the entire 
mission life time of WIRE (Bruntt & Southworth 2008). 

The four light curves from WIRE all show a double wave mod- 
ulation with a period of Prot = 4.4792 ± 0.0004 d and a peak-to- 
peak amplitude of 4mmag. We interpret this variation as being due 
to spots on the surface giving us directly the rotation period of the 
star. This is the first direct detection and was only made possible 
with the very stable photometry - on both short and long timescales 
- from WIRE. Interestingly, the amplitude and phase of the mod- 
ulation did not change significantly from 2000 to 2006 and hence 
the spot configuration seems to be extremely stable. The rotation 
period confirms an earlier indirect detection by Kurtz et al. (1994). 
Our initial assessment is that the spots need to be equatorial to ex- 
plain the double wave rotational light curve, hence - unusually for 
an Ap star - the spots are not associated with the single, visible 
magnetic pole. Further study of the rotational light variation is in 
progress. 

At much shorter periods (7 minutes) we detect the known 
dominant mode of pulsation at 2442 mHz; in addition, we have 
discovered two new frequencies located symmetrically around the 
principal mode. This triplet of frequencies is detected in the two 
runs from 2006, but not in the first two runs in 2000 and 2005. 
This is not due to differences in detection sensitivity, therefore the 
new frequencies must have appeared sometime between 2005 and 
2006. This new discovery is confirmed by our ground-based pho- 
tometry during 16 nights from SAAO, collected simultaneously 
with the two WIRE runs in 2006. The separation of the peaks is 
A/obs ~ 30.173 ± 0.004 /iHz. Based on our theoretical models of 
a Cir we interpret this as half the large separation, meaning that the 
triplet consists of modes with spherical degree of alternating parity. 
To confirm that our interpretation is correct will require that one or 
more frequencies with this separation are observed. 

The principal mode has symmetrical sidelobes with a fre- 
quency separation equal to the rotational frequency, /rot ~ 
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2.58 /iHz. We find that at the rotational phase of pulsation ampli- 
tude maximum, the pulsation phases of the three frequencies are 
equal within uncertainties. In the framework of the oblique rotator 
model this is expected when the pulsation axis and the rotation axis 
are not aligned. Thus, the three frequencies describe the change in 
observed pulsation amplitude as the observer views the deforma- 
tion of the star from different perspectives. We have used the ratio 
of the amplitudes of the rotational frequency triplet to test differ- 
ent possibilities for the spherical degree and azimuthal order of the 
principal mode. We conclude that only modes with observational 
properties similar to those of a mode with Z = 1, ?n = agree with 
the observations. This could be the dipole mode itself, or it could 
be a mode of higher degree that is distorted by the magnetic field 
such that a dipolar component dominates the observed rotational 
amplitude modulation. 

We have computed models that are consistent with the funda- 
mental parameters of a Cir using as constraints the interferometric 
radius and the effective temperature from Paper I. For the first time 
in the study of a Cir, the theoretical large separations derived from 
the models are in good agreement with the observations. The sepa- 
rations between the observed modes in the triplet /e + /i + differ 
by only 3.9 nHz and this cannot be reconciled with the theoretical 
calculations presented in this work. We have discussed that mag- 
netic effects may possibly account for this discrepancy but post- 
pone the detailed analysis to a future publication (Brandao et al., in 
preparation). 

We have compared the Hipparcos parallax of a Cir with the 
value inferred from the large separation, assuming our interpreta- 
tion above is correct. We find good agreement and this seems to 
resolve the issue of the discrepant value for a Cir as discussed by 
Matthews et al. (1999). 
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